UNCLASSIFIED 
AD  NUMBER 


DSTO  via  Itr  dtd  23  Jun  1999;  DSTO 
Library,  23  Jun  99 


THIS  PAGE  IS  UNCLASSIFIED 


I t 


E:  I € T 


i.  ^ .,  ■,  .^..  ■j.'-. 


• :/l^  .. 

■^ffii>.  '•  ■*  ■•'  *•";.*  ■■ 


|tfl<'00q-402 


teT^-ftA"0OK 


■ ' -‘  ' v!-^- 


wM0$B 


AD-C055  869 


Missile  Detection  by  Observation 
of  Differential  Motion 
in  an  Image  Sequence  (U) 


Robert  S.  Caprari 


96-00466 

mmui^ 


W-’ 


Oteders  of  the  Defence  Coi^unities^ 
Aps/ralia,  Cajma  & NZ/t^y 

hn\p  acces.^  tnis  doj2Ui^nt.  Oimri  repst' 
to  pocy^M  l^xcl^dge  (uQntryVFLS^yj 
Campl^jPark  Offices,  / / ( 

d^dfeERkA^'i  2600  mr^RAlLlA 


DTIC  CVMJTI  INSPEOTlfiB  1 


PKm 


Department. OF  Defence 


cr^Tr?XT/-tr-  a xttx  LOGY  ORGANISATION 


!('■"*•.“•  ; .V  ’V*  V J ^ V' . ' ■ . . • • 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE 
COPY  FURNISHED  TO  DTIC 
CONTAINED  A SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO 
NOT  REPRODUCE  LEGIBLY. 


Australia-Restricted 


RESTRICTED 


34. 


Copy  N<u  of  35 


Missile  Detection  by  Observation  of  Differential 
Motion  in  an  Image  Sequence  [u] 

Robert  S.  Caprari 


Land,  Space  and  Optoelectronics  Division 
Electronics  and  Surveillance  Research  Laboratory 


DSTO-RR-0018 
ABSTRACT  [U] 

This  report  investigates  a technicpic  of  passive  airborne  missile  lion 
by  analysing  a temporal  sequence  of  images  of  the  ground  in  the  region  of  the 
aircraft.  Essentially,  the  technique  identifies  as  potential  missiles,  objects  that 
move  in  a different  manner  to  the  terrestrial  background,  as  observed  in  the 
aircraft  reference  frame.  A kinernatical  analysis  is  undertaken  to  (piantify  this 
“differential  motion”.  Brief  consideration  is  given  to  the  observation  of  this  dif- 
ferential motion  by  image  processing  techniques.  The  potential  of  the  technicpie 
to  localise  the  missile  in  space  is  addressed.  The  formal  kinematics  are  applied 
to  numerical  simulations  of  guided  missile  trajectories  to  ([uantify  the  differ- 
ential motion  that  occurs,  and  determine  the  “warning  time”  between  missile 
detection  and  collision.  Several  operational  scenarios  are  investigated.  Simu- 
lation residts  are  presented,  discussed  and  appraised,  aiul  concltisions  drawn 
about  the  effectiveness  of  this  techni<itie. 
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Missile  Detection  by  Observation  of  Differential  Motion  in 

an  Image  Sequence  [in 

EXECUTIVE  SUMMARY 

This  report  addresses  the  problem  of  detection  of  a ground  launched  missile  by  an 
aircraft  without  relying  on  radar.  Although  radar  is  often  a very  effective  means  of  missile 
detection,  it  is  also  a beacon  inviting  and  facilitating  attack  by  adversaries.  One  desires  a 
method  of  observing  a missile  based  upon  light  (visible  or  infrared)  being  reflected  from  or 
emitted  by  the  missile  or  its  exhaust  gases.  An  image  of  a scene  containing  such  a missile  in 
flight  would  show  the  missile  as  a small,  bright  spot  immersed  in  a background  landscape. 
Unfortunately,  the  background  has  rapid  and  extreme  intensity  variations  which  produce 
artefacts  that  are  indistinguishable  from  mi.ssiles.  Hence,  examining  individual  images, 
even  using  the  most  sophisticated  computer  algorithms  available,  will  fail  to  distinguish 
consistently  the  real  missiles  from  false  alarms. 

One  strategy  to  overcome  this  deficiency  in  image  based  missile  detection  is  to  examine 
a serpience  of  images,  and  compare  later  images  with  earlier  ones.  Since  the  aircraft  is 
moving,  both  the  ground  an<l  juissile  images  move  in  time.  However,  because  the  missile 
is  typically  very  much  closer  to  the  aircraft,  and  also  because  the  missile  moves  relative  to 
the  ground,  the  motions  of  the  missile  and  ground  images  are  different.  The  former  effect 
is  analogous  to  the  observation  that,  when  walking  steadily  at  night,  nearby  terrestrial 
structures  (poles,  trees,  buildings)  move  into  and  out  of  view,  rviiilv  the  far  away  moon 
seems  to  follow  one  around,  because  its  position  in  the  field  of  view  doesn’t  change. 

This  report  examines  the  magnitude  of  the  motion  difference  between  the  missile  and 
ground,  as  perceived  from  the  aircraft.  It  as.sesses  whether  the  motion  is  significant  enough 
to  enable  reliable  detection  of  the  missile,  and  indeed,  whether  it  is  possible  to  estimate 
the  range  of  the  missile  from  the  aircraft  just  from  this  motion.  Essentially,  the  conclusion 
is  that  it  is  futile  to  try  to  ascertain  range  by  this  techniipie,  and  while  it  is  probably 
possible  to  detect  a missile  under  fairly  favourable  circumstances,  detection  performance 
is  not  robust  enough  to  be  considered  reliable. 
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1 Introduction 


A developing  threat  for  aircraft  operating  in  conflict  zones  is  the  proliferation  of  in- 
expensive, portable  and  lethal  infrared  guided  surface  to  air  missiles.  These  weapons  are 
certainly  accessible  to  most  conventional  armed  forces,  but  of  potentially  greater  concern 
in  a world  of  disentegrating  political  axes,  is  their  availability  to  many  rebel  armies.  There 
is  agrowing  potential  for  these  weapons  to  be  deployed  against  aircraft,  either  in  a declared 
conflict,  or  in  an  increasingly  common  peacekeeping  operation. 

Having  established  a requirement  for  defending  aircraft  against  approaching  missiles, 
one  immediately  recognises  a need  for  the  ability  to  at  least  detect  approaching  missiles, 
and  preferably  localise  the  missile  in  space.  The  most  effective  technique  of  detecting 
missiles  is  by  radar.  However,  radar  does  have  its  drawbacks.  Radio  emissions  from 
the  aircraft  provide  an  adversary  with  advance  warning  of  the  approach  of  the  aircraft, 
which  may  be  of  critical  assistance  to  the  adversary  in  mounting  an  attack  on  the  aircraft. 
Radar  is  vulnerable  to  jamming,  thus  rendering  it  inoperable  in  the  most  hazardous  cir- 
cumstances. For  reasons  of  cost,  weight,  size,  power  or  availability,  it  may  not  be  possible 
to  install  a suitable  radar  system  in  an  aircraft.  In  future  there  may  even  be  radar  seeking 
surface  to  air  missiles  in  the  arsenals  of  technologically  advanced  adversaries,  so  further 
compounding  the  problems  with  radar.  Consequently,  there  is  a demonstrable  need  for  a 
passive  missile  detection  scheme;  one  that  produces  no  electromagnetic  emissions. 

The  subject  of  this  report  is  the  proposal,  conceptual  development  and  assessment 
of  one  such  scheme  of  passive  missile  detection.  This  technique  relies  on  imaging  the 
approaching  missile  and  the  ground  below  it.  The  3-5/im  infrared  band  would  probably 
be  favoured,  because  the  missile  will  always  be  hotter  than  the  background,  either  because 
of  the  hot  plume,  or  skin  heating  by  air  friction.  Thus  the  missile  infrared  signature  will 
usually  be  conspicuous  when  it  is  immersed  in  any  realistic  background. 

Reliable  detection  of  a missile  in  individual  images  is  not  assured,  even  when  it  is 
very  bright,  because  of  the  large  intensity  variation  of  background  clutter.  Accordingly, 
the  scheme  to  be  articulated  In  this  report  exploits  the  temporal  information  available 
from  a time  sequence  of  images,  and  only  seeks  prospective  missile  candidates  that  have 
discernible  motion  relative  to  the  background.  This  “differential  motion”  is  a consequence 
of  two  factors.  One  is  that  the  missile  is  actually  moving  with  respect  to  the  background. 
The  other  is  that  the  imaging  system  is  moving  with  respect  to  both  the  ground  and 
missile,  and  the  different  distances  of  the  ground  and  missile  from  the  imaging  system 
imply  different  angular  velocities  as  seen  by  the  imager. 

This  report  addresses  the  theoretical  evaluation  of  differential  motion  from  the  specified 
motion  of  aircraft  and  missile  in  Sections  2.1  to  2.3.  Of  course,  in  practice  the  motion 
of  the  missile,  and  possibly  the  aircraft,  are  not  known.  Therefore,  differential  motion 
must  be  inferred  from  the  only  information  that  is  available,  being  the  image  sequence. 
Image  processing  operations  that  assist  in  this  endeavour  are  nominated  in  Section  2.4, 
without  elaboration.  Although  image  processing  issues  are  important  determinants  of  the 
efficacy  of  the  proposed  technique,  the  purpose  of  this  report  is  primarily  to  Investigate 
the  differential  motion  that  actually  occurs,  and  not  necessarily  its  practical  detection. 


RESTRICTED 


RESTRICTED 


DSTO-RR-0018 


Section  2.5  addresses  the  important  issue  of  the  possibility  of  utilising  the  missile  detection 
technique  that  is  developed  here  to  actually  localise  the  missile  in  space. 

Formulas  for  the  evaluation  of  differential  motion  are  applied  to  simulated  missile 
trajectories  in  Section  3.  Aspects  of  the  missile  motion  simulation  are  described,  as  is  the 
application  of  the  differential  motion  formulas  in  the  context  of  the  simulation.  Several 
simulation  results  are  presented  graphically,  and  discussed  in  terms  of  realistic  operational 
scenarios.  Section  3 concludes  with  a qualitative  and  somewhat  subjective  appraisal  of  the 
simulation  results,  based  upon  their  implications  for  useful  missile  detection  performance. 

The  report  ends  with  some  tentative  conclusions  about  the  potential  of  passive  missile 
detection  by  differential  motion  estimation  from  image  sequences. 


2 Principles  of  the  technique 

The  only  information  available  to  the  passive  missile  detection  strategy  being  inves- 
tigated is  wide  field  imagery  of  the  ground.  This  strategy  seeks  to  identify  the  preset  ■ 
of  a missile  by  discerning  its  differential  motion  with  respect  to  the  invariant  terrrestrial 
background,  as  observed  from  the  moving  reference  frame  that  is  the  aircraft. 

This  section  contains  an  exposition  of  the  kinematics  of  the  relative  motion  of  the 
ground,  aircraft  and  missile.  This  is  followed  by  analysis  that  projects  this  motion  in 
object  space  onto  the  image  plane.  The  image  plane  differential  motion  that  signifies 
the  possible  presence  of  a missile  is  then  identified.  Methods  of  observing  the  pertinent 
differential  motion  in  the  acquired  image  sequences  are  identified  in  terms  of  classical  image 
processing  operations,  to  provide  a practical  context  for  the  preceding  analysis.  Advice 
is  proffered  about  the  unsuitability  of  differential  motion  observation  for  missile  distance 
determination. 

The  following  reasonable  assumptions  simplify  the  analysis.  The  aircraft  velocity  has 
constant  direction,  and  specifically,  is  horizontal.  The  aircraft  orientation  is  invariant,  and 
specifically,  the  body  axis  of  the  aircraft  is  aligned  with  the  velocity  vector.  The  terrain 
within  the  field  of  view  is  both  flat  and  horizontal.  There  are  no  conceptual  difficulties 
associated  with  the  removal  of  any  of  these  constraints.  However,  the  intricacies  introduced 
into  the  kincmatical  analysis  by  such  a generalisation,  are  not  expected  to  invalidate  the 
insights  afforded  by  the  simple  analysis  that  is  presented  here. 


?.l  Kinematical  relationships 

Define  the  Earth  frame  O'  as  a reference  frame  that  is  fixed  on  the  Earth,  with  z'  axis 
pointing  vertically  downwards.  Define  the  aircraft  frame  O as  the  reference  frame  that  is 
parallel  to  O',  but  whose  origin  is  always  the  position  of  the  aircraft.  Note  that  under  the 
prevailing  assumptions,  the  z axis  is  always  the  optical  axis  of  the  imaging  system,  and 
its  sense  is  towards  the  field  of  the  imaging  system.  The  relationship  between  these  two 
frames  of  reference  is  illustrated  in  F'igurc  1. 
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Figure  1:  Frames  of  reference  utilised  in  the  differential  motion  analysis. 


The  positions  of  the  missile  in  the  Earth  and  aircraft  reference  frames,  respectively, 


are  (components  are  specified  in  cartesian  coordinates), 

r'm  = (*m.ym.4)  . (frame  O')  , (1) 

Tm  = (*m,ym,-»m)  , (frame  O)  . (2) 

The  positions  of  the  aircraft  in  the  Earth  and  aircraft  reference  frames,  respectively,  are 

r'a  = (®a>ya.4)  . (frame  O')  , (3) 

Po  s (0,0,0)  , (frame  O)  . (4) 


As  illustrated  in  Figure  1,  these  vectors  are  related  by  the  transformation. 


Pfn  ~ 


(5) 


or  in  cartesian  coordinates, 

Vm  = y'm-yi  y (0) 

If  frame  O is  specifically  defined  so  that  its  z axis  always  coincides  with  the  optical  axis,  and 
the  orientation  of  the  aircraft  were  such  that  an  Earth  frame  O'  could  not  be  chosen  to  be 
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aligned  with  O,  then  (5)  and  (6)  would  be  modified  by  the  inclusion  of  a general,  possibly 
time  dependent,  rotation  linear  operator,  whose  matrix  representation  is  an  orthonormal 
3x3  matrix  parametrised  by  three  suitable  Euler  angles  (eg.  roll,  pitch  and  yaw). 

Equation  (6)  represents  the  transformation  of  the  missile  and  aircraft  coordinates  in  the 
Earth  reference  frame,  into  the  position  coordinates  of  the  missile  relative  to  the  aircraft  in 
the  aircraft  reference  frame.  The  latter  frame  is  significant,  because  it  is  also  the  imaging 
system  reference  frame.  Expressing  rm  in  spherical  polar  coordinates  with 

respect  to  frame  0,  yields 

= (®m  + ym+4V^^  » 

Bm  = arctan(a!^  + € [0,7r]  , (7) 

<t>m  = arctan  • 


Assume  stigmatic,  but  in  general  distorted,  imaging  of  both  the  missile  and  the  ground. 
If  the  missile  image  is  not  resolved,  then  evidently  the  missile  image  may  be  assigned  a 
unique  position.  If  the  missile  image  is  resolved,  having  significant  spatial  extent,  one 
may  assign  the  centroid  of  the  image  as  its  unique  position.  Therefore,  one  may  sensibly 
represent  the  missile  image  position  by  a unique  point  in  both  circumstances.  The  position 
of  the  missile  image  is  completely  determined  by  only  the  angular  coordinates  of  the 
missile  itself  in  frame  O,  that  is,  Bm  and  exactly  how  will  be  elucidated  in  Section  2.2. 
Define  the  background  point  conjugate  to  the  missile,  as  that  location  on  the  ground  whose 
stigmatic  image,  in  the  absence  of  occlusion  by  the  missile,  would  coincide  with  the  missile 
image  point.  The  fact  that  formally  the  conjugate  ground  point  is  never  observed  in  the 
image  due  to  occlusion  by  the  missile,  does  not  substantively  detract  from  the  validity  of 
this  analysis,  as  will  be  explained  in  Section  2.4. 

An  immediate  corollary  of  the  discussion  of  the  previous  paragraph,  is  that  the  angular 
coordinates  in  frame  O of  the  missile  and  its  conjugate  ground  location  are  identical,  that 
is,  denoting  the  conjugate  ground  point  coordinates  by  subscripts. 


Bg  — Bjn  » 
<f>g  = <l>m  ’ 


(8) 


Conjugate  missile  and  ground  positions  are  distinguished  by  their  radial  coordinate  in 
frame  O,  that  is, 

• (9) 


The  ground  point  that  is  instantaneously  conjugate  to  the  missile  is  determined  from 
the  missile  and  aircraft  positions  by  the  following  technique.  In  the  Earth  frame  of  refer- 
ence, the  line  of  sight  from  the  aircraft  to  the  missile  has  the  parametric  equation, 

r\p)  = 4-  (r'„,  - rj,)p  , p = real  parameter  . (10) 

Resolving  (10)  into  cartesian  coordinates  yields, 

Ap)  = A + ^ 

y\p)  = ym  + (ym-yo)p»  (H) 

Ap)  = A + -^a)P  • 
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The  flat,  horizontal  ground  coincides  with  the  x'y'  plane,  so  that  in  frame  O'  it  has  the 
equation, 

«'  = 0 , (ground  plane  equation)  . (12) 

Equations  (11)  and  (12)  are  solved  simultaneously  to  reveal  the  conjugate  ground  point 
cartesian  coordinates  in  frame  O'  ^ 


In  analogy  with  (6),  the  conjugate  ground  point  cartesian  coordinates  in  frame  O 

are, 

Va  - y'g-y'a  ^ (14) 

Zg  = Zg  — z'a  , 

Expressing  (14)  in  spherical  polar  coordinates  with  respect  to  frame  O,  analogously  to  (7), 
one  obtains, 

rg  = (x^  + y^  + s^y^^  , 

0g  = arctan(x2  4- , dj€[0,)r],  (15) 

4>g  = arctan  . 


2.2  Imaging  transformation 

The  missile  detection  strategy  being  developed  here  exploits  the  displacement  in  the 
image  plane  that  arises  over  any  flnite  time  interval,  due  to  the  differential  motion  of  the 
images  of  the  missile  and  its  conjugate  ground  point.  Denote  the  stigmatic  imaging  of 
object  space  point  ro(0  to  image  plane  point  Ro(0  by  the  imaging  transformation  R 
from  a three  to  two  dimensional  space, 

Ro(f)  = R(ro(<))  , (general  imaging  transformation)  , (16) 

where  parameter  t is  the  time,  and  all  position  vectors  are  with  respect  to  the  aircraft 
reference  frame. 

The  specific  form  of  the  imaging  transformation  R is  dependent  upon  the  distortion 
of  the  optical  system.  The  case  of  perfect  imaging  (complete  absence  of  distortion)  cor- 
responds to  R being  a projective  transformation  [1,  Section  4.3].  Perfect  imaging  is  not 
practicable  for  the  extremely  wide  field  of  view  optical  systems  being  considered  for  the 
present  application,  which  ideally  have  a 180®  field  of  view.  The  problem  is  that  as  the 
field  of  view  approaches  180°,  the  size  of  the  perfect  image  increases  without  bound.  Con- 
sequently, the  periphery  of  the  field  of  view  will  transgress  the  finite  detection  region  in 
the  focal  plane,  and  will  not  be  sensed,  thereby  defeating  the  purpose  of  having  a wide 
field  optical  system. 
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A much  more  . :>'sible  imaging  transformation  would  have  sufficient  distortion  to  map 
a complete  ISO**  heid  of  vie'*/  into  a finite  region  of  the  focal  plane.  One  such  imaging 
transformation  is  that  which  actually  applies  to  some  types  offish  eye  lenses  [2][3,  p.  87,88 
and  274]  and  the  Zeiss  Pleon  lens  [4][5,  p.  220].  Define  plane  polar  coordinates  (p,  a)  in  the 
image  plane,  with  respect  the  optical-axis/focal-plane  intersection  as  origin.  The  spherical 
polar  coordinates  in  frame  O of  object  space  (r,  6,  <t>)  are  strictly  defined  with  respect  to  the 
optical-axis/entrance-pupil  intersection  ais  origin.  The  particular  imaging  transformation 
to  be  considered  is,  expressed  as  a polar  to  polar  coordinates  transformation. 

Imaging  transformation,  R : p = , 0€[0, 7r/2]  , 

a = (^  , 

where  0 is  expressed  in  radians,  and  / is  the  effective  focal  length  of  the  optical  sys- 
tem. Note  that  R depends  only  upon  the  object  space  angular  coordinates  0 and  <f>, 
as  anticipated  earlier.  Imaging  transformation  (17)  exhibits  barrel,  or  positive,  distor- 
tion. Introducing  focal  plane  cartesian  coordinates  (X,)^),  with  respect  to  the  same  origin 
as  (p,Q!),  one  has  the  final  relation. 


X = pcosa  , 
Y = psino  . 


(18) 


Combining  (7),  (17)  and  (18),  expresses  R as  a cartesian  to  cartesian  coordinates  trans- 
formation. 

Focal  length  / in  (17)  is  chosen  such  that  the  image  space,  which  is  a circle  in  the 
focal  plane,  circumscribes  the  rectangular  detector  area.  If  u)  (radians)  is  the  field  cone 
half  angle  (ie.  field  of  view  is  2a)),  PyXpx  is  the  size  of  a single  detector  element  and  tiyXrix 
is  the  size  of  the  array  of  detector  elements,  then  the  focal  length  is  chosen  as, 

/ = • (19) 


2.3  Differential  motion 

To  account  for  the  motion  of  conjugate  ground  points,  the  following  notation  is  intro- 
duced: the  ground  location  that  is  conjugate  to  the  missile  at  time  f = fi  has  the  position 
^z)  lime  t = fj.  Since  the  surface  of  the  Earth  is  essentially  rigid  (any  activity  that 
develops  in  time  is  assumed  to  occupy  a negligible  proportion  of  the  total  background),  in 
the  Earth  frame  of  reference  one  has  the  condition, 

rgih^h)  = r'g{U;t2)  Vfj  . (20) 

The  second  time  argument  in  (20)  is  thus  redundant,  and  one  may  use  the  abbreviated 
notation  in  the  Earth  reference  frame, 

r;(/,)  = r;(t,,<2)  V<2  . (21) 

In  analogy  with  (5),  one  deduces  the  relation, 

rp(fi,<2)  = r;(fi)-rUt2)  . (22) 
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By  definition,  at  an  arbitrary  time  t --  0,  the  missile  and  its  instantaneously  conjugate 
ground  point,  have  coincident  images,  that  is, 

R(r,n(0))  = R(ra(0;0))  . (23) 

For  a moving  missile,  moving  aircraft  and  stationary  background,  the  instantaneous  con- 
jugate ground  point  diverges  with  time  from  that  ground  location  that  was  conjugate  to 
the  missile  at  i = 0;  that  is,  in  general  at  time  t = T, 

t,(T!T)  . (24) 

The  differential  motion  of  interest  is  that  between  the  images  of  the  missile  and  the 
ground  location  that  was  conjugate  to  the  missile  at  the  instant  t = 0.  Denoting  the 
differential  image  displacement  that  occurs  as  a consequence  of  this  motion  during  the 
interval  from  t = 0 to  t - T by  AR(T),  the  quantity  of  interest  is, 

AR(T)  ^ R(rm(T))  - R(rg(0;T))  , (differential  motion)  . (25) 

Expressing  (25)  as  an  operation  on  position  vectors  in  the  Earth  reference  frame,  by 
substituting  (5)  and  (22),  one  obtains, 

AR(r)  = R(r;„(r)  - ri(T))  - R(ri(0)  - 4^(7))  . (26) 

Note  that  as  a consequence  of  condition  (23),  the  differential  image  displacement  has 
the  essential  property, 

AR(0)  = 0 , (27) 

that  characterises  any  vector  that  is  considered  to  be  a displacement  that  develops  over 
time. 


2.4  Image  processing  interpretation 

The  interpretation  of  differential  image  displacement  (25)  in  terms  of  image  processing 
operations  is  clarified  by  invoking  (23),  and  expanding  (25)  to  become, 

AR(T,  = [R(rm(T))  - R(r,,(0))]  - [R(r,(0;  T))  - R(r,(0;  0))]  . (28) 

The  quantity  [R(rm(T’))  - R(rm(0))]  is  measured  by  identifying  the  missile  image 
point  in  each  of  two  images  corresponding  to  the  instants  i = 0 and  t = T,  and  evaluating 
the  vector  displacement  between  the  two  points.  This  is  essentially  a manifestation  of  the 
classical  target  identification  problem  in  image  processing.  If  there  is  more  than  one  missile 
present,  or  if  there  are  additional  erroneous  missile  identifications,  one  must  determine 
the  association  between  individual  targets  in  the  initial  and  final  images.  Associativity  is 
allocated  according  to  one  of  the  established  methods  of  statistical  decision  theory. 

In  principle,  the  quantity  [R(rp(0;r))-R(rj(0‘,0))]  may  be  measured  in  an  analogous 
manner,  where  now  it  is  the  conjugate  ground  point  at  f = 0 that  is  to  be  identified 


RESTRICTED 


7 


RESTRICTED 


DSTO-RR-0018 


in  the  two  images.  In  practice  this  is  infeasible,  because  this  point  is  occluded  by  the 
missile  in  the  original  image,  so  that  one  has  no  reference  with  which  to  identify  this 
point  in  the  final  image.  However,  the  fact  that  the  background  is  extended,  and  that 
the  variation  with  position  of  the  background  displacement  between  initial  and  final  im- 
ages is  slow  compared  with  the  characteristic  length  of  the  background  image  structure 
(that  is,  the  image  distortion  due  to  optical  system  motion  is  slight  compared  with  the 
intrinsic  variation  in  the  background  image),  offers  a feasible  method  of  measuring  the 
displacement  [R(rj(0;T))-  R(rj(0;0))]. 

One  may  select  a region  of  background  surrounding  the  missile  image  in  the  initial 
image  to  be  a “template".  Assuming  that  the  size  of  the  missile  image  is  small  compared 
with  the  image  distortion  due  to  optical  system  motion,  one  selects  the  size  of  the  back- 
ground region  to  be  much  greater  than  the  missile  image  size.  One  then  optimally  matches 
the  final  image  to  this  template,  by  a correlation  operation.  The  resulting  displacement 
of  the  correlation  maximum  is  interpreted  as  being  the  best  estimate  of  the  displace- 
ment [R(rj(0;T))  — R(rj(0;0))].  Although  the  slight  distortion  between  tlie  two  images, 
and  the  different  positions  of  the  small  missile  images  relative  to  the  background  in  the  two 
images,  will  degrade  the  equivalence  between  the  correlation  maximum  displacement  and 
the  displacement  being  sought,  the  stated  assumptions  ensure  that  this  perturbation  is 
marginal.  This  procedure  is  essentially  a manifestation  of  the  classical  image  registration 
problem  in  image  processing.  Since  image  registration  is  inherently  a nonlocal  operation, 
it  is  not  materially  affected  by  the  occlusion  of  conjugate  ground  points  by  the  missile,  as 
noted  earlier. 

In  summary,  the  classical  image  processing  operations  of  target  identification  and  image 
registration,  possibly  supplemented  by  the  decision  theoretic  operation  of  associativity, 
measure  the  quantities  that  determine  the  differential  displacement  AR(T)  via  (28).  This 
establishes  another  essential  property  of  the  definition  of  AR(T)  (25);  the  existence  of 
an  identifiable  procedure  for  its  measurement  from  the  only  available  information,  being 
sequences  of  images. 


2.5  Suitability  of  differential  motion  for  missile 
localisation 

Although  the  detection  of  a missile  is  itself  a commendable  achievement,  it  is  desirable 
also  to  determine  the  missile  displacement  from  the  aircraft.  Since  differential  motion 
observation  mandates  the  use  of  an  imaging  sensor,  the  direction  of  the  missile  in  the 
aircraft  frame  is  dcducible.  The  angular  coordinates  of  the  missile  with  respect  to  the 
imager  optical  axis  are  determined  from  the  missile  image  plane  coordinates  by  the  inverse 
of  the  imaging  transformation  deflned  by  (16).  Specifically,  one  would  apply  the  inverses 
of  (18)  and  (17)  in  that  order,  to  convert  the  missile  image  coordinates  (Am.Vm)  to  the 
object  space  angular  coordinates  {0mi<t>m)‘ 

The  only  scope  for  inferring  missile  range  is  from  the  observed  differential  motion.  It 
is  apparent  that  the  differential  motion  between  the  missile  and  conjugate  ground  point 
images  depends  upon  the  missile  trajectory,  and  indeed,  the  ground  topography.  In  the 
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absence  of  knowledge  of  these  two  factors,  especially  the  former,  it  is  not  possible  to  make 
a credible  determination  of  missile  range  from  differential  motion  observation. 

For  instance,  there  exists  a family  of  stealthy  trajectories  along  which  a hypothetical 
missile  could  approach  the  aircraft  without  the  conjugate  ground  point  ever  being  per- 
turbed. Hence,  there  would  never  be  any  differential  motion  to  observe,  and  the  missile 
would  never  be  detected.  More  realistic  is  the  example  of  a missile  whose  trajectory  in  the 
aircraft  frame  is  always  directly  towards  the  aircraft.  This  is  actually  a good  approxima- 
tion for  the  trajectory  arising  from  a proportional  navigation  guided  missile  approaching 
an  aircraft  flying  with  constant  velocity.  In  such  a situation,  the  position  of  the  missile 
image  is  invariant,  and  any  differential  motion  that  is  observed  is  entirely  a consequence 
of  the  motion  of  the  ground  relative  to  the  aircraft.  In  both  of  these  circumstances,  the 
first  rather  artificial,  but  the  second  quite  realistic,  it  is  not  even  in  principle  possible  to 
infer  the  missile  range  by  observing  its  differential  motion,  since  the  missile  range  docs 
not  influence  the  differential  motion. 

For  a general  missile  trajectory,  the  missile  range  does  influence  the  differential  motion 
in  a trajectory  dependent  manner.  If  the  missile  trajectory  and  ground  topography  were 
both  known,  then  there  would  be  a formal  solution  for  the  missile  range  in  terms  of 
the  differential  displacement  magnitude.  However,  in  the  absence  of  independent  sensing 
facilities,  the  missile  trajectory  in  particular  is  not  known,  so  it  is  not  possible  to  reliably 
determine  missile  range  from  differential  motion  observation  in  realistic  circumstances. 

In  summary,  the  use  of  an  imaging  sensor  renders  the  missile  angular  coordinates 
deducible  once  the  missile  has  been  detected  by  observing  differential  motion  in  the  image 
sequence.  In  contrast,  the  missile  range  can  not  be  reliably  inferred  from  differential 
motion  observation  in  the  variable  circumstances  encountered  in  practice. 


3 Simulation  of  missile  flight  and  detection 

Having  established  the  principles  of  missile  detection  by  observation  of  the  differential 
displacement  AR(T),  it  remains  to  assess  the  effectiveness  of  the  technique  in  scenarios 
that  will  be  encountered  In  practice.  This  is  accomplished  by  simulating  the  trajectory  of  a 
guided  missile  from  launch  to  engagement  of  the  aircraft,  while  simultaneously  calculating 
the  resulting  differential  motion. 

In  this  section,  the  missile  trajectory  simulation  is  described  superficially,  and  with 
the  assistance  of  the  exposition  of  Section  2,  a description  of  the  computation  of  differ- 
ential motion  is  presented.  Results  of  particular  simulations  are  presented,  discussed  and 
appraised. 


3.1  Features  of  the  missile  trajectory  simulation 

Realistic  modelling  of  powered  and  guided  missile  flight  in  the  atmosphere  is  a formid- 
able undertaking.  The  present  investigation  was  expedited  by  utilising  a missile  simulation 
program  furnished  by  other  researchers  [6,  7].  This  simulation  program  was  modified  as 
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outlined  later  in  this  section,  and  incorporated  as  an  essential  element  in  a differential 
motion  calculation  program. 

The  missile  simulation  is  nominally  a 3 degrees-of-freedom  approximation  of  the  mis- 
sile dynamics,  in  which  the  three  missile  position  coordinates  are  the  only  free  kinematical 
parameters.  To  account  for  the  thrust  direction,  the  missile  orientation  is  constrained 
to  be  tangential  to  the  trajectory.  This  condition  implies  two  nonholonomic  constraints 
(constraints  expressed  as  differential  equations;  in  this  case  the  orientation  coordinates 
determined  by  the  time  derivatives  of  the  centroid  coordinates  [8,  Section  1.3]),  since  ori- 
entation is  specified  by  two  angles  (elevation  and  azimuth).  A third  holonomic  constraint 
(constraint  expressed  as  an  algebraic  equation)  is  that  the  roll  angle  of  the  missile  is  in- 
variant. These  constraints  account  for  the  three  “missing”  degrees  of  freedom  compared 
with  the  actual  situation.  However,  to  account  for  the  drag  on  the  missile,  the  inclina- 
tion of  the  missile  from  the  trajectory  tangent  is  allowed  to  be  a free  pair  of  parameters 
(elevation  and  azimuth).  Effectively,  the  simulation  is  a “3  and  a bit”  degrees-of-freedom 
approximation  of  the  missile  trajectory. 

Ironically,  the  presence  of  nonholonomic  constraints  in  the  simplified  model  of  missile 
dynamics  implies  that,  even  in  principle,  the  equations  of  motion  probably  do  not  have 
a formal  solution.  In  the  absence  of  formal  analytic  solutions,  the  present  simulation 
numerically  solves  the  equations  of  motion  by  finite  difference  methods. 

The  approximate  nature  of  the  missile  model  inevitably  causes  a significant  discrep- 
ancy between  the  simulated  trajectory  and  the  actual  trajectory.  The  severity  of  this 
discrepancy  has  not  been  quantified,  but  the  author  believes  that  this  simulation  is  at 
least  qualitatively  correct,  and  that  the  quantitative  errors  are  not  so  extravagant  as  to 
render  the  simulation  unrepresentative  of  realistic  missile  motion. 

The  missile  characteristics  utilised  in  the  present  simulations  have  been  selectively 
modified  from  those  existing  in  the  original  simulation.  Variation  of  mass  and  thrust  with 
time,  external  dimensions  for  evaluation  of  drag,  fuel  load  for  range  determination  and 
the  guidance  law  (proportional  navigation  in  this  case),  are  all  representative  of  an  actual 
missile.  However,  unlike  any  real  missile,  in  these  simulations  there  is  no  maximum  limit 
on  the  angular  velocities  or  angular  accelerations  of  either  the  missile  body  or  the  seeker 
head;  effectively  the  missile  and  seeker  have  zero  rotational  inertia.  Furthermore,  there 
is  also  no  maximum  limit  on  the  angular  displacement  of  the  seeker  axis  from  the  missile 
axis. 

Although  the  original  missile  simulation  algorithm  did  impose  bounds  on  the  achievable 
missile  dynamics,  these  bounds  were  deliberately  removed  from  the  missile  model  used  for 
the  present  simulations  to  enable  a more  extensive  test  of  the  detection  algorithm.  It 
transpires  that  with  the  manoeuvrability  bounds  set  at  realistic  values,  there  is  typically 
a large  region  on  the  ground  within  which  a launched  missile  fails  to  engage  the  aircraft, 
according  to  the  original  missile  simulation.  By  removing  the  selected  constraints,  the 
engagement  launch  envelope  has  its  “voids”  significantly,  but  not  entirely,  filled  in.  Since 
one  is  most  interested  in  assessing  the  performance  of  the  missile  detection  algorithm, 
and  not  the  missile’s  ability  to  engage  its  target,  it  is  beneficial  to  manipulate  the  missile 
characteristics  to  maximise  the  opportunities  of  testing  the  detection  algorithm  (such  a 
test  requires  the  missile  to  intercept  the  aircraft). 
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Removal  of  the  missile  dynamics  bounds  almost  certainly  increases  the  inaccuracy  of 
the  missile  trajectory  simulation  beyond  that  which  already  afflicted  the  original  simula- 
tion. No  attempt  was  made  to  gauge  the  severity  of  degradation  of  accuracy.  The  author 
is  of  the  opinion  that  this  further  compromise  in  accuracy  does  not  materially  affect  the 
qualitative  assessment  of  the  results  of  these  investigations.  However,  the  existence  of  this 
further  source  of  Inaccuracy  does  reinforce  the  need  for  caution  in  the  interpretation  of 
the  results.  A prudent  and  justifiable  policy  to  adopt,  would  seem  to  be  only  to  draw 
qualitative  conclusions  from  one’s  assessment  of  the  results. 


3.2  Computation  of  differential  motion 

Having  solved  the  missile  m-:*tion,  knowing  the  uniform  aircraft  motion  and  knowing 
that  the  ground  is  rigid,  sufficient  information  is  now  available  to  compute  the  differential 
motion  AR(T)  from  the  kinematical  considerations  enunciated  in  Section  2. 

The  missile  trajectory  computation  furnishes  cartesian  coordinates  of  the  missile  po- 
sition vector  r'^  in  the  Earth  frame.  The  cartesian  coordinates  of  the  aircraft  position 
vector  Ta  in  the  Earth  frame  are  defined  by  the  motion  of  the  aircraft.  Having  chosen  a 
suitable  time  interval  T for  the  displacements  to  develop,  the  differential  n.otion  at  arbi- 
trary time  t = 0 is  evaluated  by  the  direct  application  of  (26).  As  explained  in  Section  2.2, 
imaging  transformation  R is  expressed  as  a cartesian  (R^)  to  cartesian  (R^)  coordinates 
operation  by  combining  (7),  (17)  and  (18). 


3.3  Method  of  presentation  of  simulation  results 

Aspects  of  the  computation  of  the  differential  motion  corresponding  to  a single  missile 
launch  position  have  been  addressed  in  Sections  3.1  and  3.2.  To  obtain  full  coverage  of  the 
ground,  these  simulations  are  repeated  on  a relatively  fine  grid  (40  to  100  metres  spacing) 
of  launch  sites  on  the  ground.  The  grid  of  launch  sites  extends  beyond  the  range  of  the 
missile  to  encompass  the  complete  threat  region. 

The  results  of  this  substantial  simulation  are  best  interpreted  if  one  evaluates  a single 
critical  parameter  as  a scalar  field,  with  the  argument  being  the  launch  position.  Undoubt- 
edly, the  most  critical  parameter  is  the  warning  time,  defined  as  the  time  that  elapses 
between  the  differential  motion  algorithm  first  deciding  that  a missile  is  present,  and  the 
missile  actually  engaging  the  aircraft.  Results  of  the  current  simulations  are  displayed  as 
quantised  intensity  plots  of  the  warning  time  scalar  field  (often  referenced  as  warning  time 
“maps”  in  this  document). 

All  simulations  correspond  to  a constant  velocity  aircraft.  The  simulation  program  ac- 
tually does  accomodate  climbing  (positive  elevation  from  horizontal)  or  descending  (neg- 
ative elevation)  aircraft,  despite  the  kinematics  established  in  Section  2.1  assuming  hori- 
zontal flight.  Axes  of  the  parallel  reference  frames  O and  O'  are  oriented  with  respect  to 
the  aircraft  as  illustrated  in  Figure  1,  that  is;  a;,  x' — from  right  to  left;  y,  y' — from  rear  to 
front;  z,  z' — downwards  along  the  optical  axis.  The  aircraft  velocity  is  along  the  positive 
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y'  direction.  A non-zero  elevation  is  achieved  by  rotating  the  aircraft  body  and  its  velocity 
about  the  x axis,  in  a clockwise  sense  through  the  prescribed  elevation  angle. 

Regardless  of  the  aircraft  elevation,  there  is  reflection  symmetry  of  the  kinematics  in 
the  plane.  Accordingly,  the  warning  time  field  is  only  computed  and  displayed  for 
launch  positions  with  x*  > 0;  the  field  for  z'  < 0 is  implied  by  refiection  in  the  y'  axis. 
In  all  of  the  displayed  maps,  the  aircraft  position  at  the  instant  of  missile  launch  is  such 
that  z'  = 0,  y'  = 0 and  z'  < 0.  The  position  coordinates  of  the  maps  are  the  z'  and  y' 
coordinates  of  the  missile  launch  site,  and  the  launch  position  specification  is  completed 
by  setting  z'  = - 1 metres. 

The  algorithm  decides  that  it  has  detected  a missile  according  to  the  following  simplis- 
tic, but  reasonable,  heuristic  decision  process.  The  differential  displacement  vector  AR(T) 
is  expressed  in  units  of  horizontal  and  vertical  detector  element  dimensions,  and  its  Eu- 
clidean norm  is  compared  with  a prescribed  threshold.  If  this  threshold  is  exceeded  for 
a prescribed  number  of  consecutive  images  in  the  sequence,  the  algorithm  decides  that  it 
has  detected  a missile.  The  time  at  which  the  missile  is  detected  is  noted,  but  the  missile 
trajectory  computation  continues  until  engagement  of  the  aircraft,  this  time  also  being 
noted.  The  warning  time  is  defined  as  the  time  that  elapses  between  the  missile  being 
detected  and  the  missile  engaging  the  aircraft. 

The  displayed  maps  indicate  several  intervals  of  warning  time  by  different  intensities. 
Black  shading  generally  corresponds  to  zero  warning  time,  that  is,  the  algorithm  fails  to 
detect  the  missile  at  any  time  before  engagement.  The  one  possible  exception  to  this  rule 
is  the  region  in  the  vicinity  of  the  origin  of  the  z'y'  plane.  The  missile  trajectory  simulation 
program  is  unstable  (ie.  the  computed  motion  diverges  rapidly  from  a realistic  trajectory) 
in  this  region,  and  one  adopts  the  conservative  approach  of  assigning  a zero  warning  time 
within  this  region.  The  extent  of  this  region  is  chosen  to  be  a circle  that  subtends  an  angle 
of  10**  from  the  vertical  at  the  aircraft  initial  position.  White  shading  always  indicates 
failure  of  the  missile  to  engage  the  aircraft,  regardless  of  whether  the  missile  was  detected 
before  the  trajectory  simulation  was  terminated.  Since  all  maps  extend  beyond  the  range 
of  the  missile,  they  all  have  white  regions  on  their  periphery.  Regions  of  lighter  shading  on 
individuai  maps  signify  longer  warning  time  intervals.  Except  for  black  and  white  regions, 
there  is  in  general  no  warning  time  equivalence  between  regions  of  the  same  shading  in 
different  maps.  Accordingly,  one  must  exercise  caution  when  comparing  different  maps;  in 
particular,  one  should  not  necessarily  attach  the  same  significance  to  regions  of  identical 
shading  in  different  maps. 

The  elapsed  time  interval  T for  the  development  of  displacements  from  the  differential 
motion,  is  always  chosen  to  be  an  integral  multiple  of  the  time  interval  separating  con- 
secutive images  in  the  sequence  of  images.  Several  aircraft  dispositions  are  considered  in 
the  simulations.  For  each  disposition,  all  imaging  and  algorithm  parameters  are  invariant, 
except  for  T.  Warning  time  maps  are  computed  for  three  different  values  of  T for  every 
aircraft  disposition.  This  facilitates  assessment  of  the  crucial  influence  that  T exerts  on 
the  attainable  warning  times.  There  is  comparative  flexibilty  in  the  choice  of  T,  since  it 
is  essentially  a software  parameter,  and  is  not  influenced  by  the  optical  system  apparatus. 
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3.4  Discussion  of  warning  time  maps 

Computed  warning  time  maps  corresponding  to  seven  realistic  operational  scenarios 
are  displayed  in  Figures  2 to  8 on  pages  18  to  45.  Features  of  the  warning  time  maps 
are  identified  and  discussed  in  this  section,  without  consideration  at  this  stage  of  the 
implications  of  the  simulation  results.  The  latter  undertaking  is  reported  in  Section  3.5. 

The  warning  time  maps  of  Figure  2 are  representative  of  a high  speed  aircraft  in  transit. 
Only  the  time  interval  T within  which  the  differential  displacement  develops  differs  among 
the  three  maps.  A significant  void  in  the  warning  time  map,  within  which  the  missile  fails 
to  engage  the  aircraft,  is  evident  as  the  irregularly  shaped  white  region  intruding  into  an 
otherwise  convex  shaded  region  of  successful  engagement.  Such  voids  are  present  to  varying 
degrees  in  other  warning  time  maps.  The  extent  to  which  the  warning  time  maps  as  a 
whole,  and  the  voids  in  particular,  accurately  reflect  real  missile  performance,  is  dependent 
upon  the  quality  of  the  missile  trajectory  calculation.  As  indicated  in  Section  3.1,  this 
computation  is  not,  nor  is  it  designed  ^o  be,  particularly  accurate.  Comparing  the  three 
maps  demonstrates  that  within  much  of  the  region  of  successful  engagement,  the  warning 
time  increases  noticeably  with  increasing  integration  time  T,  as  expected,  and  consistent 
with  most  other  scena'ios  investigated. 

The  warning  time  maps  of  Figure  3 are  representative  of  very  fast  flight  at  moderate 
altitude.  In  this  situation,  the  void  in  the  warning  time  map  particularly  severely  curtails 
the  engagement  envelope.  To  some  extent  this  may  be  due  to  lack  of  anticipation  in 
the  choice  of  initial  missile  direction;  the  missile  is  always  launched  directly  towards  the 
aircraft,  regardless  of  the  aircraft  speed.  Once  again,  detection  of  the  missile  improves 
dramatically  with  integration  time  T. 

The  warning  time  maps  of  Figure  4 are  indicative  of  a transport  aeroplane  having  just 
taken-off  (note  the  positive  elevation).  The  upward  elevation  of  the  aircraft,  together  with 
its  low  altitude,  are  such  that  a missile  launched  from  far  behind  the  aircraft  will  approach 
the  aircraft  nearly  horizontally,  so  the  missile  will  be  beyond  the  imaging  system  field  of 
view.  Consequently,  such  a missile  will  never  be  detected,  regardless  of  the  integration 
time  T.  This  inference  is  confirmed  by  the  three  warning  time  maps,  all  of  which  exhibit 
V lishing  warning  time  for  launches  that  are  far  enough  behind  the  aircraft. 

The  warning  time  maps  of  Figure  5 are  illustrative  of  an  aircraft  on  a jungle  reconnais- 
sance mission.  Within  almost  all  of  the  engagement  launch  envelope  there  is  no  warning. 
Further  simulation  demonstrated  that  if  the  field  of  view  is  increased  to  180®,  then  the 
region  of  no  warning  is  reduced  by  a moderate  amount,  but  still  dominates  the  warning 
time  map.  These  warning  time  maps  are  anomalous  in  the  sense  that  the  warning  time 
does  not  increase  monotonically  with  integration  time  (map  (c)  has  the  longest  integra- 
tion time,  but  the  shortest  warning  time).  This  result  cautions  one  against  assuming  that 
missile  detection  performance  necessarily  increases  with  integration  time  with  no  limit, 
even  for  very  long  integration  times.  It  doesn’t! 

The  warning  time  maps  of  Figure  6 represent  very  low  altitude,  high  speed  flight.  For 
these  maps  the  optical  system  field  of  view  has  had  to  be  increased  to  a full  hemisphere 
to  obtain  any  detection  at  all.  The  fact  that  a slight  increase  in  field  of  view  from  160° 
to  180°  has  transformed  the  missile  detection  performance  from  complete  failure  to  quite 
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satisfactory,  attests  to  the  critical  influence  of  the  optical  system  field  of  view  on  missile 
detection  by  very  low  altitude  aircraft. 

The  warning  time  maps  of  Figure  7 are  indicative  of  a transport  aircraft  during  long 
distance  transit.  The  performance  of  the  differential  motion  missile  detection  technique 
shows  substantial  improvement  with  integration  time. 

The  warning  time  maps  of  Figure  8 are  illustrative  of  a helicopter  cruising  horizontally. 
This  simulation  considers  the  helicopter  to  be  oriented  horizontally,  instead  of  having 
its  nose  inclined  downwards  to  achieve  forward  motion,  as  for  actual  helicopter  flight. 
These  warning  time  maps  are  all  distinguished  by  the  complete  absence  of  voids  in  the 
engagement  launch  envelope.  The  warning  time  that  is  achieved  by  the  differential  motion 
observation  increases  with  integration  time. 


3.5  Appraisal  of  warning  time  maps 

One  seeks  to  assess  the  performance  of  the  differential  motion  missile  detection  tech- 
nique on  the  basis  of  the  warning  time  maps  of  Figures  2 to  8.  Apart  from  the  rather 
exceptional  scenario  represented  by  Figure  5,  where  the  missile  was  usually  at  the  periph- 
ery of  the  field  of  view  and  the  integration  times  were  relatively  long,  the  warning  time 
increases  with  integration  time.  Therefore,  each  scenario  will  be  assessed  on  the  basis  of 
its  longest  T warning  time  map,  being  map  (c)  in  each  figure.  One  will  arbitrarily,  but 
realistically,  presume  that  warning  times  of  less  than  2 seconds  are  totally  inadequate, 
those  of  2-4  seconds  are  marginally  useful,  and  those  above  4 seconds  are  satisfactory. 

Figure  2(c)  exhibits  a reasonable  area  of  successful  detection,  but  there  remains  a 
significant  region  of  unsuccessful  detection  in  the  forward  direction.  The  performance  of 
the  missile  detection  algorithm  would  probably  improve  demonstrably  with  an  even  longer 
integration  time.  However,  in  the  absence  of  confirmation  of  this  hypothesis,  one  concludes 
that  for  this  scenario,  missile  detection  fails  for  an  unacceptably  large  region. 

Figure  3(c)  indicates  quite  successful  missile  detection  within  most  of  the  engagement 
launch  enveope.  The  only  exception  is  a small  region  towards  the  exteremity  of  the 
engagement  envelope  in  the  forward  direction.  Increasing  the  integration  time  could  render 
the  missile  detection  successful  essentially  everywhere  within  the  launch  envelope  for  this 
scenario. 

Figure  4(c)  is  characterised  by  the  majority  of  the  region  within  the  launch  envelope 
corresponding  to  unsuccessful  missile  detection.  Furthermore,  for  much  of  the  region  be- 
hind the  aircraft,  the  failure  to  detect  the  missile  is  probably  due  to  the  missile  flying 
above  the  extremity  of  the  field  of  view.  Increasing  the  field  of  view  to  a full  hemisphere 
would  probably  improve  the  missile  detection  performance  in  the  backward  direction  sig- 
nificantly. Further  enhancement  of  the  missile  detection  capability  would  accrue  from 
increasing  the  integration  time.  However,  on  the  basis  of  existing  evidence,  missile  detec- 
tion by  differential  motion  is  unsatisfactory  for  this  scenario. 

All  throe  warning  time  maps  of  Figure  5 emphatically  confirm  that  missile  detection 
by  diffcrontial  motion  is  most  unsatisfactory  in  this  scenario.  These  computations  suggest 
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that  there  are  no  advantages  to  be  gained  by  increasing  the  integration  time  even  further. 
The  limiting  factor  in  the  performance  of  the  missile  detection  algorithm  over  much  of  the 
region  is  the  field  of  view.  Expanding  the  field  of  view  to  180®  yields  a moderate  increase 
in  warning  time  over  much  of  the  region,  however,  the  overall  performance  of  the  missile 
detection  strategy  remains  inadequate  for  this  scenario. 

Figure  6(c)  shows  satisfactory  missile  detection  performance  within  most  of  the  en- 
gagement launch  envelope,  an  outcome  that  is  only  enabled  by  using  a 180®  field  of  view. 
The  regions  for  which  missile  detection  is  inadequate  would  probably  be  reduced  in  size 
by  an  increase  in  integration  time.  Hence,  missile  detection  by  differential  motion  is  qtiitc 
successful  for  this  scenario. 

Figure  7(c)  exhibits  successful  missile  detection  within  much  of  the  engagement  enve- 
lope, although  there  remains  an  intolerably  large  region  in  the  forward  direction  within 
which  missile  detection  is  unsuccessful.  Increasing  integration  time  should  have  a mod- 
erately beneficial  impact  on  missile  detection  performance,  but  on  the  basis  of  only  the 
presented  evidence,  the  performance  is  inadequate. 

Figure  8(c)  exhibits  approximately  equal  areas  within  the  engagement  envelope  corre- 
sponding to  each  of  the  inadequate,  marginal  and  satisfactory  missile  detection  categories. 
One  could  expect  a modest  increase  in  warning  times  with  an  increasing  integration  time, 
but  even  with  this  enhancement,  the  overall  missile  detection  performance  would  still  be 
inadequate  for  this  scenario. 


4 Conclusion 


The  differential  motion  missile  detection  simulations  that  have  been  presented  in  this 
document,  while  not  exhaustive,  have  elucidated  some  important  factors  that  influence 
the  performance  of  missile  detection.  Maximising  the  field  of  view  has  considerable  bene- 
fits, which  are  especially  dramatic  for  low  flying  aircraft.  Increasing  the  integration  time 
within  which  the  differential  displacement  is  allowed  to  develop  is  the  other  factor  that 
dramatically  improves  performance,  under  almost  all  circumstances. 

Nevertheless,  on  the  basis  of  the  present  calculations,  one  would  conclude  that  differ- 
ential motion  observation  does  not  provide  one  with  a reliable  missile  detection  strategy 
over  a diversity  of  scenarios.  The  number  of  simulations  that  have  been  conducted  is  in- 
sufficient to  convey  an  indication  of  any  pattern  in  the  success  of  the  algorithm  in  different 
scenarios,  if  indeed  such  a pattern  even  exists.  Given  the  nature  of  the  differential  motion 
technique,  its  missile  detection  performance  depends  critically  on  the  missile  trajectory. 
Therefore,  it  would  probably  be  futile  to  seek  such  a hypothetical  pattern,  since  it  would 
almost  certainly  vary  with  different  missile  trajectories,  thus  not  being  of  general  validity. 

There  certainly  is  scope  for  more  detailed  investigation  of  the  effect  of  various  factors 
on  missile  detection  by  differential  motion  observation.  Given  the  reservations  expressed 
about  the  validity  of  the  missile  trajectory  computation  employed  in  the  present  inves- 
tigations, perserverance  with  the  current  simulations  is  not  recommended  for  discerning 
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subtle  patterns  in  the  outcome,  which  may  be  destroyed  by,  or  even  arise  from,  simula- 
tion inaccuracy.  However,  the  author  is  of  the  opinion  that  the  qualitative  assessment 
of  the  technique  stated  here  will  not  be  overturned  by  further  accurate  simulation:  that 
differential  motion  observation  has  the  potential  to  detect  approaching  missiles  only  under 
specific  circumstances  that  are  too  restrictive  for  it  to  be  considered  a candidate  for  the 
primary  missile  detection  strategy. 

Regardless  of  whether  differential  motion  observation  is  deemed  a satisfactory  means 
of  detecting  missiles,  it  has  been  argued  in  Section  2.5  that  it  is  certainly  not  a reliable 
method  of  determining  missile  range.  This  detracts  further  from  the  attractiveness  of 
differential  motion  observation  as  the  primary  missile  detection  strategy.  It  seems  that  at 
best,  observation  of  differential  motion  should  only  be  used  to  augment  another  generally 
superior  missile  detection  strategy,  if  indeed  there  is  a need  for  such  a supportive  role. 

A further  comment  on  the  scope  of  these  investigations  may  be  appropriate.  The 
advice  proffered  above  relates  only  to  the  possibility  of  ranging  or  detection  of  missiles  by 
observation  of  differential  motion.  In  no  way  should  it  be  construed  as  having  relevance 
to  any  other  passive  ranging  and  detection  techniques,  other  infrared  techniques,  or  other 
image  processing  techniques  that  seek  to  achieve  the  same  capability. 
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Figure  2 

(next  three  pages) 

List  of  common  parameter  values. 


Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  fim 
Pixel  period  in  the  y direction  = 50.0  ftm 
Field  of  view  cone  half  angle  = 80.0* 

Frame  rate  = 25.0  s“* 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 983.0  m 
Upwards  elevation  = 0.0* 

Aircraft  ground  speed  = 1075.0  km/h 


Warning  tintes  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-3 

• 3-6 

• 6-9 

• 9-12 

• 00  (no  engagement;  white  region) 
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MISSILE  DETECTION  MAP 
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Figure  2:  (a)  Time  interval  for  differential  displacement  to  develop  = 2 frames. 
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Figure  2:  (b)  Time  interval  for  differential  displacement  to  develop  = 6 frames. 
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Figure  2:  (c)  Time  interval  for  differential  displacement  to  develop  = 10  frames. 
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(next  three  pages) 

List  of  common  parameter  values. 


Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  /im 
Pixel  period  in  the  y direction  = 50.0  /im 
Field  of  view  cone  half  angle  = 80.0® 

Frame  rate  = 25.0  s“* 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 1764.0  m 
Upwards  elevation  = 0.0® 

Aircraft  ground  speed  = 1356.0  km/h 


Warning  times  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-2 

• 2-4 

• 4-6 

• 6-8 
• 8-10 

• 00  (no  engagement;  white  region) 
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Figure  S:  (a)  Time  interval  for  differential  displacement  to  develop  = 3 frames. 


RESTRICTED 


RESTRICTED 


map  05 

DSTO-RR-0018 


MISSILE  DETECTION  MAP 


10000 


transverse  launch  position  (metres) 


Figure  3:  (b)  Time  interval  for  differential  displacement  to  develop  = 8 frames. 
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Figure  S:  (c)  Time  interval  for  differential  displacement  to  develop  = 16  frames. 
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Figure  4 

(next  three  pages) 

List  of  common  parameter  values. 

Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  /xm 
Pixel  period  in  the  y direction  = 50.0  /xm 
Field  of  view  cone  half  angle  = 80.0° 

Frame  rate  = 25.0  s“’ 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 397.0  m 
Upwards  elevation  = 8.0° 

Aircraft  ground  speed  = 173.0  km/h 

Warning  times  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-2 

• 2-4 

• 4-6 

• 6-8 
• 8-10 

• 00  (no  engagement;  white  region) 
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Figui'e  4-'  (b)  Time  interval  for  differential  displacement  to  develop  = IS  frames. 
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Figure  4';  (c)  Time  interval  for  differential  displacement  to  develop  = 18  frames. 
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Figure  5 

(next  three  pages) 

List  of  common  parameter  values. 

Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  ^m 
Pixel  period  in  the  y direction  = 50.0  fim 
Field  of  view  cone  half  angle  = 80.0® 

Frame  rate  = 25.0  s”* 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 277.0  m 
Upwards  elevation  = 0.0® 

Aircraft  ground  speed  = 159.0  km/h 

Warning  times  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-1 
• 1-2 

• 2-3 

• 3-4 

• 00  (no  engagement;  white  region) 
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Figure  5:  (a)  Time  interval  for  differential  displacement  to  develop  = 12  frames. 
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Figure  5:  (b)  Time  interval  for  differential  displacement  to  develop  = 25  frames^ 
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Figure  5:  (c)  Time  interval  for  differential  displacement  to  develop  = 50  framec. 
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List  of  common  parameter  values. 


Number  of  pixels  in  the  x direction  =:  512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  fim 
Pixel  period  in  the  y direction  = 50.0  (tm 

Field  of  view  cone  half  angle  = 90.0** 
Frame  rate  = 25.0  s“' 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  s:  168.0  m 
Upwards  elevation  ==  0.0** 

Aircraft  ground  speed  = 580.0  km/h 


Warning  times  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-3 

• 3-6 

• 6-9 

• 9-12 

• 12-15 

• 00  (no  engagement;  white  region) 
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Figure  6;  (a)  Time  interval  for  differential  displacement  to  develop  = 5 frames. 
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Figure  6:  (h)  Time  interval  for  differential  displacement  to  develop  = 10  frames. 
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Figure  6:  (c)  Time  interval  for  differential  displacement  to  develop  = 20  frames. 
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List  of  common  parameter  values. 


Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  fim 
Pixel  period  in  the  y direction  = 50.0  /im 
Field  of  view  cone  half  angle  = 80.0** 

Frame  rate  = 25.0  8“* 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 2307.0  m 
Upwards  elevation  = 0.0® 

Aircraft  ground  speed  = 322.0  km/h 


Warning  times  (seconds)  corresponding  to  intensity  levels: 

• 0 (no  warning;  black  region) 

• 0-3 

• 3-6 

• 6-9 

• 9-12 


• 00  (no  engagement;  white  region) 
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Figure  7;  (a)  Time  interval  for  differential  displacement  to  develop  = S frames. 
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Figure  7:  (b)  Time  interval  for  differential  displacement  to  develop  = 7 frames. 
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Figure  7:  (c)  Time  interval  for  differential  displacement  to  develop  — 12  frames. 
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List  of  common  parameter  values. 

Number  of  pixels  in  the  x direction  = 512 
Number  of  pixels  in  the  y direction  = 512 
Pixel  period  in  the  x direction  = 50.0  /im 
Pixel  period  in  the  y direction  = 50.0  /im 
Field  of  view  cone  half  angle  = 80.0® 

Frame  rate  = 25.0  s“* 

Detection  threshold  displacement  = 2.00  pixels 
Required  consecutive  missile  detections  = 4 
Initial  height  of  aircraft  = 824.0  m 
Upwards  elevation  = 0.0® 

Aircraft  ground  speed  = 119.0  km/h 

Warning  times  (seconds)  corresponding  to  intensity  levels; 

• 0 (no  warning;  black  region) 

• 0-2 

• 2-4 

• 4-6 

• 00  (no  engagement;  white  region) 
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Figure  8:  (a)  Time  interval  for  differential  displacement  to  develop  = 5 frames. 
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Figure  8:  (b)  Time  interval  for  differential  displacement  to  develop  = 15  frames. 
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Figure  S:  (c)  Time  interval  for  differential  displacement  to  develop  = 25  frames. 
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